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ABSTRACT 

Strong Balmer absorption lines and the lack of Ha and [Oil] emission lines signify that E+A 
galaxies are post-starburst systems. Recent studies suggest that E+As may undergo the transi- 
tion from the 'blue cloud' to the 'red sequence' and eventually migrate to red sequence early- 
type galaxies. An observational validation of this scenario is to identify the intervening galaxy 
population between E+As and the red-sequence. Motivated by recent findings with Galaxy 
Evolution Explorer (GALEX) that an unexpectedly large fraction of early-type galaxies ex- 
hibit UV-excess (i.e. blue UV - optical colours) as a sign of recent star formation (RSF), we 
investigate the possible connection of the UV-excess galaxies to E+As. In particular, we ex- 
amine the Fundamental Plane (FP) scaling relations of the currently largest sample of ~ 1,000 
E+As selected from the SDSS and ^^20,000 morphologically-selected SDSS early-type galax- 
ies with GALEX UV data. The FP parameters, combined with stellar population indicators, 
reveal a certain group of UV-excess early-types that bridges between E+As and quiescent red 
galaxies. The newly identified galaxies are the post-starburst systems characterized by UV- 
excess but no Ha emission. This is essentially a conceptual generalisation of "E+A", in that 
the Balmer absorption line in the "E+A" definition is replaced with UV - optical colours that 
are far more sensitive to RSF than the Balmer lines. We refer to these UV-excess galaxies as 
"E+a" galaxies (named after "E+A"), which stands for elliptical ("E") galaxies with a minor- 
ity of A-type ("a") young stars. The species are either (1) galaxies that experienced starbursts 
weaker than those observed in E+As (1 ~ 10 per cent of E+As, "mild E+As") or (2) the 
products of passively evolved E+As after quenching star formation quite a while ago (~ 1 
Gyr, "old E+As"). We suggest that the latter type of E+a galaxies (i.e. old "E+As") represents 
the most recent arrival to the red sequence in the final phase of the E+A to red early-type 
transition. 

Key words: galaxies: fundamental parameters - galaxies: elliptical and lenticular, cD - galax- 
ies: formation - galaxies: evolution - galaxies: starburst - ultraviolet: galaxies 



1 INTRODUCTION 

The study of galax y evolution began to flourish after 
iDressler & GunnI Jl983h discovered galaxies with unusual spectra 
that showed characteristics of both elliptical galaxies and A-type 
stars, now referred to as "E+A" galaxies. E+As are considered as 
post-starburst systems because of strong Balmer absorption lines, 
sign of recent starbursts within ~ 1 Gyr, and a lack of [Oil] and 
Ha emission lines inclicating a suddenly quenched star forma tion 
(e.g. lDressler & Gunn|[l983l : iPoggianti et al]ll999l : lGotoll2004 . 

How star formation in E+As has started and quenched is a 
topic that has been studied by many. One plausible explanation is 
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that E+As result from some physical mechanis ms in cluster en- 
vironments, such as ram-pressure stripping ( e.g. ISpitzer & Baadd 
ll95ll : lGunn & Gotl'l972';'Fuiita & G otol2004h . galaxy-cluster tidal 
interaction (Bryd & Valtonen 1990), or galaxy harassment (e.g. 
iMoore et al.lll996 ). However, the fact that a significant fraction of 
E+As are found in field environments implies that cluster-related 
physical m echanisms are not enough to explain the origin of a ll 
E+As (e.g. IZabludoff et al.lll996l : iBlake et al.ll2004l: lGotdl20o"5bh . 
Another possible explanation for the E+A phenomena is starburst 
galax ies with very thick optical depth caused by dust obscura- 
tion dCouch & Sharpies! 1 19871 : ISmail et al.|[T99^ : IPoggianti & wiil 
I2000I) . but radio and IR observations have detected little star for- 
mation obscured by dust in E+As (Galazl |2000| : iMiller & Ower] 
l2001. : .Goto..2004.) . Recently. iGotg C2005b.) showed that E+As have 



2 Y. Choi, T. Goto and S.-J. Yoon 



closer companions than other galaxies do, providing evidence that 
a dynamical merger / interaction could be the physical origin of 
field E+As. This scenario is also supported by the dynamically dis- 
turbed morphologies of E+As tTran et al . 2004; Yang et al. 2 0041 : 
lYamauchi & Gotd l2005l : iLiu et alj 1200% . While there are some 
plausible scenarios, the origin of E+As is still not well understood. 

E+As are the key populations for understanding the galaxy 
evolution because they are believed to be in a transition phase from 
a 'blue cloud' to a 'red sequence'. Base d on galaxy proper ties such 
as colour profiles and scaling relations, lYang et al] j2008l) suggest 
that E+As are the remnants of galaxy-galaxy interactions/mergers 
and ultimately evolve into early-type galaxies. The typical migra- 
tion time from the gas-rich 'blue cloud' to the gas-poor 'red se- 
quence' is about 1.5 Gyr according to the UV properties of E+As 
jKavirai et ai]|2007ah . UV light is useful for tracking E+A migra- 
tion because it is very sensitive to the young stellar population. The 
recent st ar formation ( RSF) signature in UV light continues about 
1.5 Gyr dYTetal .120051) . According to recent GALEX observations, 
a significant fraction of early-type galaxies in the local Universe 
have enhanced UV flux (i.e. blue UV - optical colours) as a sign 
of young stellar population (e.g. Yoon et al. 2004a; Yi et alj|2005l ; 
ISalim et a l. 2005; Donas et al. 2007; Kaviraj et al. 2007b). The rel- 
atively weak RSF galaxies might have arisen mainly through two 
channels: (a) the residual star formation in gas-poor galaxies, or (b) 
the fading away of a young stellar population in galaxies that un- 
derwent a violent starburst. UV-excess early-types appear to be the 
most recent arrival at the red sequence in the final phase of the blue 
cloud to red sequence transition, and may thus be linked to E+As. 

Early-type galaxies can be described by three observational 
parameters: their effective radius (re), effective mean surface 
brightness (^) within re, and central velocity dispersion (cr). These 
parameters are uni fied in a two-dimensional manifold, the Funda - 
mental Plane (FP) jPiorgovski & Davisll987l;lDressler et al.ll987h . 
Even though the FP has intrinsic scatter ( |j0rgensen et aljl99d) . the 
thin FP of early-type galaxies implies that there are well defined 
scaling relations between their global parameters. The FP can be 
expressed such that log = a log a + f3 fj. + y, where i?e is 
the effective radius in kpc. Supposed that early-type galaxies are 
virialized homologous systems with a constant mass-to-light ratio, 
the virial theorem predicts a = 2.0 and (3 = 0.4. However, the ac- 
tual, measured coefficients are different from the theoretical expec- 
tations (called the "tilt" of the FP). The FP is a useful tool for ex- 
plori ng evolutionary connections among various classes of g alaxies 
(e.g. lGenzel et al.ll200ll : lTacconi et alJl2002l ; lLee et"Z|[2007t) . The 
purpose of this paper is to investigate the evolutionary path from 
E+As to quiescent galaxies via UV-excess galaxies by comparing 
their FPs. 

The paper is organized as follows. A description of the sample 
used in this study is given in § 2. We examine the FPs of our sam- 
ples and discuss the implication of our results in § 3. We will report 
an evolutionary connection between E+As and UV-excess galaxies. 
Finally, § 4 gives the summary. 

Except where stated otherwise, we adopt the ACDM cosmol- 
ogy with {h, n™, f^A) = (0.7, 0.3, 0.7). 
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Figure 1. SDSS pipeline velocity dispersions vs. our measurements for 
673,668 galaxies in the SDSS DR6 catalog. The two dashed vertical lines 
ai'e the lower (70 km s~^) and the upper (430 km s~^) limits of the ve- 
locity dispersion in our sample. The instrumental dispersion of the SDSS 
spectrograph sets the lower limit (69 km s^^). The upper limit is given 
above which none of template spectra fits the observed ones for most of the 
galaxies. Red circles are the 377 E+As in common, satisfying 70 km s~^ 
< (T < 430 km s^^ in both measurements. 



as lGotd JioOTh . i.e. H(5 EW > 4 A, [Oil] EW > -2.5 A, Hq EW 
> -3.0 A, and a signal-to-noise ratio > 10 pixeP^ in the r-band. 
Since velocity dispersions were measured for only 415 E+As in the 
SDSS pipeline, we independently measure the velocity dispersions 
of 673,668 galaxies in the SDSS DR6. Our method is similar to that 
utilized by the SDSS spectroscopic pipeline; we first mask out pos- 
sible emission line regions, and then fit the rest of the spectra with 
a combination of eigen spectra with varying velocity dispersions. 

Figure[T]compares our velocity dispersion measurements with 
those from the SDSS pipeline, showing reasonable one-to-one 
matches. Red circles are the 377 E+As in common, satisfying 70 
km s~^ ^ cr ^ 430 km s~^ in both measurements. We restrict our 
sample to galaxies with 70 km s~^ ^ a 430 km s~^ because the 
instrumental dispersion of the SDSS spectrograph is 69 km s^^ per 
pixel, and none of template spectra matches the observed ones for 
most of the galaxies with a of 430 km s^^ or higher. Dashed verti- 
cal lines mark the lower (70 km s^^) and upper (430 km s^^) lim- 
its on the velocity dispersion of our sample. Among the 1,284 E+A 
candidates from the SDSS DR6, we will use 1,021 E+As which 
satisfy the velocity dispersion criterion. This is the largest and the 
most homogeneous E+A sample ever investigated. 



2 SAMPLE 
2.1 E+A galaxies 

This study uses the largest sample of 1,284 E+As selected from 
the SDSS DR6. The selection criteria of the E+As are the same 



2.2 GALEX / SDSS early-type galaxies 

The GALEX / SDSS early-type galaxy sample is constructed as fol- 
lows. First, we make use of ~ 180,000 morpholog i cally- selected 
SDSS DR4plus early-type galaxies from lChoi et al.l ( l2007h. In this 
morphological classification method first presented by Park & Choil 
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Figure 2. Distributions of 2-band absolute mag, A/^, and redsliift, z, of tlie 
sample used in this study. Solid lines are for E+A galaxies and dotted lines 
are for the early-type galaxies with GALEX UV information. 



( |2005|) . u — r colour, radial gradient in Q — i colour, and co ncen- 
tration are the criteria for classification. IPark & Choil l l2005l) esti- 
mated the completeness and reliability of this classification method 
reaching about 91 per cent for SDSS galaxies brighter than r-band 
Petrosian mag, r^et < 15.9 and about 88 per cent at Vpet < 17.5. 
Then, we match the SDSS early-type galaxies with UV detections 
in the 2,241 GALEX fields from the GALEX DR4, imaged in the 
Medium Imaging Survey (MIS), the Deep Imaging Survey (DIS), 
and Nearby Galaxy Survey (NGS) modes. The cross-matching is 
performed with an angular matching radius of 3 arcsec. For the 
multiple GALEX / SDSS matches within the given matching radius 
(~ 5 per cent of all matches), we select the UV source with the 
smallest angular distance from the SDSS galaxy. The final catalog 
consists of 23,539 early-type galaxies with near-UV (NUV) infor- 
mation (~ 15 per cent of the entire early-type galaxy sample in 
Choi et al. 2007). 

To give basic information of our sample, Figure |2] provides 
the distributions of 2-band absolute mag, Mz, and redshift of E-l-A 
and early-type galaxies. The early-type sample shows very simi- 
lar Mz distribution to that of E-l-A, suggesting their similar stellar 
mass ranges. However, due to the GALEX UV detection limit, the 
smaller galaxies are selectively missing in the higher redshift bins, 
resulting in early-types being confined to lower-redshift regime 
compared to E-l-A. In order to avoid the selection bias, the com- 
parative analysis between early-type galaxies and E-l-As will be re- 
stricted to the early-type sample with z < O.I. Also note that our 
volume-limited early-type sample will have redshift ranging from 
0.025 to 0.09. 

Figure [3] illustrates the power of UV detecting young stel- 
lar populations. Most of the quiescent and UV-excess early-type 
galaxies are so narrowly distributed in g — r colour (~ 0.2 mag) 
that UV-excess early-type galaxies cannot be distinguished from 
the quiescent ones. In contrast, early-type galaxies with RSF are 
clearly separated from the quiescent ones in NUV — r colour. 
Note that the FUV flux is often affected by the UV upturn while the 
NUV flux is less sensitive to th e phenomenon dYoon et al.ll2004bl : 
lLeeetalJ|200^; iRee etalj|20OT") . Thus, the NUV flux is a bet- 
ter tracer for RSF activities in early-type galaxies. This is why 
NUV — r colour is used as a selection criterion for UV-excess 
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Figure 3. The H5 vs. g — r diagram (upper panel) and the H<5 vs. NUV — r 
diagram (lower panel) for QST, RSF I, and RSF 11. Both RSF I and RSF II 
ai'e indistinguishable from QST in the g — r colour distribution. Contrary to 
g—r colour, there is a separation between RSF I /II and QST in the NUV — 
r colour distribution. A gap between QST and RSF \/Ilm NUV — r colour 
is an artifact of their different selection criteria. 



early-type galaxies. The criterion for UV-excess early-ty pe galax- 
ies is rather conservatively defined as NUV — r < 5.5 jYi et al.l 
l2005l : lKavirai et al.ll2007tt ISchawinski et al.ll2007h . The threshold 
comes from the NUV — r colours of nearby giant early-type galax- 
ies with the strongest UV upturn but no sign of RSF a ctivities (e.g 
iBurstein et al]|l988l: luee et aI.ll2005L iRee et al.ll2007l) . Therefore, 
we can not rule out the possibility that this threshold excludes some 
of the less massive RSF galaxies, but, a better threshold than this 
should await more precise information about the dependence of UV 
upturn properties on galaxy mass. 

The selected UV-excess early-type galaxies are further divided 
into the following two types by using NUV ~ r colour and another 
prime RSF indicator. Ha emission line: (a) Ongoing weak star for- 
mation galaxies (RSF I) that are characterized by the presence of 
Hq emission line with an RSF signature in NUV — r colour, and 
(b) Post-weak-star-formation or bygone weak star formation galax- 
ies (RSF II) that show a sign of RSF in NUV — r colour but no 
Ha emission. The rest of the galaxies are UV-dead quiescent galax- 
ies that show no evidence of new-born stars, and are referred to as 
"QST" . 

The |Pai-k& Choil ( l2005h morphological classification method 
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Figure 4. Distributions of QST (top), RSF I (middle), and RSF II (bottom) 
galaxy subsamples in the u — r colour vs. radi al g — i colour gra dient 
plane. This is the same as Fig. 1 (lower panel) of IPark & Choil llOOSi) , but 
with early-type galaxies having GALEX UV information. The exact shape 
of each line is as a function of galaxy rp^t as denoted in the bottom panel. 



Table 1. Four different RSF modes in early-type galaxies. Early-type galax- 
ies in our sample are classified by using Ha emission line and NUV — r 
colour 



RSF mode 


Ha 
Emission line 


NUV - r 
colour 


Class 


Quiescent mode 


No 


Red 


QST 




(> 0.5 A) 


(6.5 < , ^ 7.5) 




Ongoing weak SF 


Weak 


Blue 


RSF I 




(-9A<,^-3A) 


«5.5) 




Bygone weak SF 


No 


Blue 


RSF II 




(> 0.5 A) 


« 5.5) 


"E+a" 


Ongoing starburst 


Strong 




SHq 




(«; -9 A) 







2.0 



1 .5 



1 .0 



0.5 



Orthogonal fit 



(5.42) 



(5.54) 



(6.02) 



(7.69) 



(14.75) 




Direct fit 



r-band (30.43) 
9-band\ (33.90) 
I I I I 



QST RSF I RSF II SHa E + A 



performs well in general. However, since the strategy relies on 
galaxy colours, the UV-excess early-type galaxies may be selec- 
tively missed in this scheme. It is important to specify how this 
affect our results. F igure [4] is the same as Fig. 1 (lower panel) of 
IPark & Choil ilOO^ but with early-type galaxies having GALEX 
UV information. It is clear that QST and RSF II galaxies are fairly 
secure in this scheme, while about a quarter of RSF I galaxies ap- 
pear to be more consistent with the late-type morphology. This 
agrees with the result based on visual inspections in Figures [T2lll3l 
andll4l Since our results are based mainly on the QST and RSF II 
galaxies, the colour-based classification scheme does not affect our 
conclusion. 

In addition, we select 1,548 early-type galaxies with a strong 
Hq emission line from the SDSS DR6 with the following criteria: 

• Inverse concentration index, rgo/rso < 0.33 

• b/a (axis ratio) in the r-band ^ 0.6 

• Hof emission line EW ^ — 9 A 

• S/N in the r-band > 10 

» AGN-free on the BPT diagram jSaldwin et aP Il98ll : 
iKewlev et alj200ll : lGotdl2005 j) 

• Exclusion of 36 disc bulges by eye inspection 

We refer to these galaxies as "SHa" (strong Ha) and they do not 
necessarily have the GALEX UV data. The selection criteria for 



Figure 5. The evolution corrected coefficients of all galaxy types obtained 
from the direct fitting (circles) and the orthogonal fitting (squares) methods. 
Galaxies in different RSF modes have different FP slopes (a). The solid and 
dashed lines are the coefficients in the r-band and the 3-band, respectively. 
The coefficients show wavelength dependency in all types of galaxies. Note 
that the FP coefficients vary with fitting methods. The direct fit gives the 
steepest slope for the QST FP and the shallowest for the SHa FP in both 
bands. The orthogonal fit yields the steepest slope for the SHa FP in the r- 
band and the RSF II FP in the g-band. The figures in parentheses represent 
the dihedral angles between our g-band planes and the virial plane in units 
of degree. 



the different types of galaxies introduced in this section are sum- 
marized in Table[T] 



3 THE FUNDAMENTAL PLANES OF E+A GALAXIES 
AND GALEX /SDSS EARLY-TYPE GALAXIES 

3.1 Construction of the Fundamental Planes 

The three observed parameters for the FP analysis are obtained as 
follows. We set the effective radius to = {h/af^'^Tdev, where 
b/a IS the ratio of the minor and major axes and rdeu is an ef- 
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Table 2. Coefficients of tlie FPs of tlie five types of early-type galaxies in the r-band and the g-band. Errors are estimated from the bootstrap method (10,000 
resamplings). The rms of the direct fit and the rms of the orthogonal fit are the scatter around the plane in the log fi e direction and in the perpendicular 
direction, respectively. We apply the evolution coiTection following the evolution parameter Q of lBemardi et alj<2003cl) . 



r-band 



-band 



Direct Fit 



QST 


1.25 ±0.01 


0.3 10 ±0.003 


-8.33 ±0.07 


0.058 


1.20 ±0.01 


0.307 ± 0.003 


-8.46 ± 0.08 


0.047 


RSFl 


1.13 ±0.03 


0.254 ± 0.006 


-6.91 ±0.16 


0.166 


1.08 ±0.03 


0.238 ± 0.006 


-6.68 ±0.17 


0.166 


RSFII 


1.18±0.02 


0.254 ± 0.005 


-7.06 ±0.11 


0.160 


1.12±0.02 


0.239 ± 0.005 


-6.83 ±0.11 


0.156 


SHa 


0.59 ±0.03 


0.200 ±0.008 


^.60 ±0.19 


0.191 


0.57 ±0.03 


0.222 ± 0.006 


-5.14±0.15 


0.194 


E-hA 


0.73 ±0.03 


0.206 ± 0.004 


^.98 ±0.12 


0.184 


0.68 ± 0.02 


0.209 ± 0.003 


-5.05 ±0.10 


0.170 


Orthogonal Fit 


















QST 


1.61 ±0.02 


0.306 ±0.004 


-9.09 ±0.08 


0.053 


1.60 ±0.02 


0.301 ±0.004 


-9.23 ±0.09 


0.055 


RSFl 


1.72 ±0.05 


0.286 ±0.007 


-8.82 ±0.22 


0.072 


1.75 ±0.05 


0.266 ± 0.007 


-8.71 ±0.23 


0.078 


RSF 11 


1.75 ±0.03 


0.270 ± 0.005 


-8.62 ±0.1 3 


0.075 


1.75 ±0.03 


0.249 ± 0.005 


-8.43 ±0.14 


0.080 


SHa 


1.80±0.12 


0.230 ±0.011 


-7.83 ±0.34 


0.129 


1.76±0.12 


0.251 ±0.009 


-8.30 ±0.33 


0.130 


E-hA 


1.20 ±0.04 


0.222 ± 0.005 


-6.35 ±0.14 


0.108 


1.05 ±0.04 


0.217 ±0.003 


-6.06 ±0.12 


0.104 


Evolution-Conected Direct Fit 


















QST 


1.19±0.01 


0.311 ±0.003 


-8.23 ±0.07 


0.014 


1.12±0.01 


0.308 ± 0.003 


-8.32 ±0.07 


0.035 


RSFl 


1.10±0.03 


0.258 ±0.006 


-6.92 ±0.15 


0.159 


1.03 ±0.03 


0.244 ± 0.006 


-6.71 ±0.16 


0.156 


RSF 11 


1.14±0.02 


0.258 ±0.005 


-7.06 ±0.11 


0.150 


1.06 ±0.02 


0.244 ± 0.004 


-6.82 ±0.11 


0.142 


SHq 


0.58 ±0.03 


0.211 ±0.008 


^.79 ±0.18 


0.183 


0.55 ± 0.03 


0.230 ±0.006 


-5.27 ±0.14 


0.182 


E+A 


0.70 ±0.03 


0.210 ±0.004 


-5.02 ±0.12 


0.175 


0.63 ±0.02 


0.211 ±0.003 


-5.03 ±0.09 


0.161 


Evolution-Con'ected Orthogonal Fit 


















QST 


1.51 ±0.02 


0.306 ± 0.003 


-8.87 ±0.07 


0.051 


1.46 ±0.02 


0.300 ± 0.003 


-8.94 ± 0.07 


0.053 


RSFl 


1.62 ±0.05 


0.284 ± 0.007 


-8.60 ± 0.20 


0.071 


1.62 ±0.05 


0.264 ± 0.007 


-8.44 ±0.21 


0.076 


RSFII 


1.65 ±0.03 


0.270 ±0.005 


-8.43 ±0.12 


0.074 


1.63 ±0.03 


0.249 ± 0.005 


-8.20 ±0.13 


0.079 


SHa 


1.68±0.12 


0.232 ±0.011 


-7.63 ±0.31 


0.127 


1.59±0.11 


0.251 ±0.008 


-7.95 ± 0.30 


0.128 


E-hA 


1.12±0.04 


0.223 ± 0.004 


-6.20 ±0.13 


0.105 


0.95 ±0.03 


0.217 ±0.003 


-5.87 ±0.11 


0.102 



fective radius from the de Vaucouleurs fit jde Vaucouleurs|[T948h . 
We do this in orde r to convert an ellipti cal aperture to an effec- 
tive circular radius tBemardi et al.l2003al) . The effective mean sur- 
face brightness within re is given by = rudev + 2.5 log(27rre) 
— K(2:) — 10 log(2:-l-l), where mdeu i s the appar ent magnitude ob- 
tained by the de Vaucouleurs' r^^'* fit dde Vaucouleurs 1948) and z 
is the redshift. We take into account the K-correction as described 
bv lBIanton et alj (12003 ) and the cosmological (l-l-z)* dimming ef- 
fect (e.g. lTolmaniri93C ). We correct the velocity dispersion for the 
aperture effect of the SDSS using (Teo r = (-^^7lr^)'^'°'^g-eat, where 



r fiber is 1.5 arcsec and is in arcsec ( |j0rgensen et al J 19951) . 

We determine the FP coefficients (a, /3, 7) of each type 
of galaxies (QST, RSF I /II, SHa, and E-l-A) by using the two 
schemes; the least square direct fitting method and the least square 
orthogonal fitting method. By performing the robust fit, we reduce 
the effect of outliers. Direct fitting yields the coefficients that min- 
imize the thickness of the plane in the log Re direction, while the 
orthogonal fitting algorithm minimizes the summation of residuals 
perpendicular to the plane. The plane coefficients are obtained from 
simple algebraic calculations. For the orthogonal fitting, we use a 
following covariance matrix 

for which the eigenvector corresponding to the smallest eigenvalue 
is the normal vector of the plane we are seeking. The smallest 
eigenvalue represents the rm s scatter of the plane (e.g. lSaglia et al.l 
I2OOII : iBemardi et alj[20033) . For the direct fitting method, on the 



other hand, we use the same algebra as lBemardietai] ( l2003cl) . The 
uncertainties in the coefficients are estimated via 10,000 bootstrap 
resamplings of the data. The best-fit coefficients and their errors are 
presented in Table[2] 

For a sanity check, we make a combined subset of QST and 
week-SF (i.e. RSF I and RSF II) galaxies and compare the FP co- 
efficients to those of Bernardi et al. (2003a). We assume that by 
putting SHq galaxies aside the subset mimics their sample. The de- 
rived coefficients from both the direct fit a nd the orthogonal fit ar e 
about 5 ~ 10 per cen t larger than those of lBemardi et alj ( l2003d ). 
iBemardi et al.l (l2003b) showed that the corrections for the evolu- 
tion and selection effects tend to reduce the coefficient a. They de- 
termined the evolution parameter Q using the maximum likelihood 
method in the g-band (0.85) and in the r-band (1.15). After tak- 
ing weak passive luminosity evolution effect into account, we find 
the c oefficients for t he com bined sample in very good agreement 
with iBemardi et al] ( l2003cl) . For example, the evolution corrected 
coefficients a of our combined sample are 1.07 for the direct fit 
and 1.50 for the orthogonal fit in the p-band. The Q-correction has 
been applied to all classes of our sample (QST, RSF I /II, SHa, and 
E-l-A). The slopes of these subclasses bec ome shallower and con- 
sistent with the coefficients suggested by iBernardi et al.l ( l2003cl) 
within the coefficient measurement errors. The evolution corrected 
coefficients are listed in Table|2] Note that we do not apply the cor- 
rection for the selection effect which is no more straightforward 
due to complications arising from the cross matching between two 
datasets (SDSS and GALEX). Nevertheless, the bias in a due to 
this effect is relatively ins ignificant compared to the evolution ef- 
fect dBemardi et al.ll2003cl) and hence does not affect our results. 



6 Y. Choi, T. Goto and S.-J. Yoon 





Figure 6. Edge-on view of the r-band FP for five different types of early-type galaxies. Red circles are E+A (every panel) and black circles are QST (upper-left 
panel), RSR I (upper-right), RSF II (lower-left), and SHa (lower-right), respectively. Lines represent the best fit of each sample in the edge-on projection of 
the virial plane. The distribution of E+A is most analogous to that of strong Ha galaxies. RSF II also exhibits relatively strong similarities to E+A. 



Figure [5] shows the evolution corrected coefficients a of all 
galaxy types obtained from the direct fit (circles) and the orthogo- 
nal fit (squares). Th e well-known wav e length dependency of the 
FP coefficients (e.g. Pahre et al. 1998; Scodeggio et al. I [1991) IS 
weak (e.g Bemardi et al. 2003c; Colless et al. 2001) but apparent. 
It has been pointed out that the choice of data fitting method affects 
the FP coefficients (e .g. iSagha et al.ll200ll ; iBemardi et Zll2003cl ; 
iD'Onofrio eT^bOOSl) . The systematic difference among various 
fittin g methods is most likely caused by the scatter of the FP (e.g. 
iLa B arbera et al. 2000). Given that there is no standard method for 
deriving the coefficients, it is important to choose the method that 
describes the given data most properly. The coefficients a are es- 
pecially important values when they are used as indicators of the 
FP tilt. To test which fitting methods work as more reliable and 
robuster indicators of the FP tilt, we examine a dihedral angle be- 
tween our plane and the virial plane. The dihedral angles in degree 
are denoted in parentheses (Fig.|5}. While the coefficients a from 
the orthogonal fit do not show a systematic tendency, those from 
the direct fit behave well as a function of the dihedral angles, in 
the sense that the coefficients a decrease as the dihedral angles in- 
crease. For the orthogonal fit, the rank of a among different galaxy 
types can even change depending on the band in use. For instance. 



the orthogonal fit gives the steepest slope (i.e. the greatest a) for 
SHq in the r-band and for RSF II in the g-bwA, respectively. It is 
also worth noting that while the orthogonal fit has the advantage 
of treating the FP parameters symmetrically, the method yields co- 
efficients with larger uncertainties, especially for SHa. Taken to- 
gether, the coefficients a from the direct fit more successfully de- 
scribe the relative features of the FP slopes among different galaxy 
types. Hence, we use the coefficients a based on the direct fitting 
method in the following analysis. 



3.2 Comparison of the Fundamental Planes of E+As and 
Early-type Galaxies in Different RSF Modes 

lYang et alj ( l2008b measured the coefficients of the FP of 16 E+As 
in the Gunn r and the Johnson B using the orthogonal fitting 
method. Their best-fit FP has a = 1.13 ± 0.10 and /3 = 0.24 ± 0.07 
for the Gunn r, and a = 1.09 ±0.08 and /3 = 0.23 ±0.06 for 
the Johnson B. The largest E+A sample (~ 60 times larger than 
that used in Yang et al. 2008) allows us to estimate the E+A FP 
with higher accuracy. Our E+A FP based on the orthogonal fitting 
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Figure 7. Same as Fig.|6] but in the g-band. Since the g-band is more sensitive to the young stellar population than the r-band, the separation among different 
types becomes more obvious. 



method is given by 

log_Re = 1.20(1.12) log cr + 0.222(0.223) M - 6.35(6.20) , 
logi?e = 1.05(0.95) log cr + 0.217(0.217) M - 6.06(5.87) 

in the r- band and g-band, respectively. The evolution corrected co- 
efficients appear in parentheses. When corrected for the evolution 
effect, a becomes 1.12 in the r-band and 0.95 in the g-band, while 
P rema ins nearly constan t in both bands. Although the bandpass 
used in I Yang et al.l feOOSl) differs slightly from the SDSS r-band 
and g-band, our E+A FP obtained from the orthogonal fit agrees 
well with that o f I Yang etalJ dlOOS') within the uncertainties regard- 
less of whether the evolution effect is corrected or not. The E+A FP 
calculated by the direct fitting method is 

log_Re = 0.73(0.70) log cr + 0.206(0.210) ^ - 4.98(5.02) , 
logi?e = 0.68(0.63) log cr + 0.209(0.211)^ - 5.05(5.03) 

in the r-band and the g-band, respectively. 

Figures |6] and |7] show the distribution of E-l-A on the edge-on 
view of the virial plane, as compared to four other types of galaxies; 
QST, RSF I, RSF II, and SHa. The differences between E+A and 
the other three types of galaxies are more clearly seen in the g- 
band (see Fig.|7](, which is more sensitive to the presence of young 



stars than the r-band. The comparison shows that E+A lies on a 
plane that is different from those of QST, RSF I, and RSF II in 
terms of slope, position, and scatter. As mentioned in Section [JTI 
the coefficients a obtained from the direct fitting method are used 
when comparing the FP slopes. It is evident that E+A is absent in 
the upper right part and more extended toward the lower left part of 
the edge-on projection. Among the four types of galaxies, the QST 
FP is the most distinct from the E+A FP, both in slope and in scatter. 
The FP slope of E+A corresponds to ~ 60 per cent of the QST 
slope in both the r-band and the g-band, and the rms scatter (i.e. the 
thickness) of the E+A FP is nearly three times greater than that of 
the QST FP. In terms of position, the E+A FP extends more toward 
the lower left part as compared to the QST FP. Comparing RSF I to 
E+A, we find a similar slope, but a significantly different position. 
RSF I is distributed over a narrow range, whereas E+A occupies a 
relatively long and broad region on the edge-on projection of the 
virial plane. On the contrary, RSF II is similar to E+A in both the 
slope and position. The SHa FP shows a stronger similarity to that 
of E+A with slightly larger (1.15 times) scatter. The scatter around 
the SHa FP is the largest, and the scatter of the E+A, RSF I / II, and 
QST follow it in the order. 
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Figure 8. The normalized distributions of log Re, log a, and 11 of E+A and QST in the g-band. The black and red histograms present the distributions of 
the FP parameters of QST and E+A, respectively. Solid curves are fitted Gaussian functions. In all parameters, E+A shows a larger standard deviation and a 
smaller mean than QST. The lower right panel shows the Kormendy relation of E+A and QST with 2.2 ^ log cr < 2.4. The Kormendy relation of E+A is 
steeper than that of QST, leading to the tilt of the E+A FP. 



3.3 Origin of the tilt and shift of the E+A Fundamental Plane 

In Figure [S] we compare the distributions of three FP parameters 
of QST (black) and E+A (red) in order to investigate the physi- 
cal cause of the shift and tilt of the E+A FP. Each distribution is 
normalized by its maximum value. The solid curves are the Gaus- 
sian fits to the histograms, the standard deviation and the mean of 
which are summarized in Table [3] Compared to QST, E+A has 
smaller means and larger standard deviations in all parameters of 
interest. The smaller effective radius and more luminous surface 
brightness (i.e. smaller ^ val ue) of E+A are likely due to centrally 
concentrated star formation ( ICaldwelletaTlll999l) . A Starburst in 
the galaxy centre increases the central luminosity, leading to both 
a decreased effective radius and increased mean surface brightness 
within the effective radius (left panels of Fig.[8}. As a consequence, 
E+A migrates toward the lower left corner of the virial plane pro- 
jection, as compared to QST (see Fig. |6] and Fig.|7}. The smaller 
velocity dispersion of E+A is also involved in the shift of the E+A 
FP. According to the downsizing scenario, less massive galaxies 
have experienced more active star formation in the recent epoch, 
and so one would expect E+A to be less massive and thus have a 
smaller velocity dispersion than QST. 

In order to test whether the smaller velocity dispersion of E+A 
is an intrinsic property, we construct a volume-limited sample of 
0.025 ^ 2 ^ 0.09 and ^ -20.5. The criteria include galaxies 
with a wide range of luminosities and a maximum number of E+As 
in the sample. Figure |9] (left panel) shows that the volume-limited 
sample exhibits the same trends as seen in the entire E+A sample, 
namely, a smaller mean and lager standard deviation in the log cr 



Table 3. The means and the standard deviations of the distributions for log 
Re, log (T, /X, and log cr + 0.2/x of all types in the g-band. 
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log 


(J 






log (T + 0.2/i 




mean 


std 


mean 


std 


mean 


std 


mean std 


QST 


0.740 


0.274 


2.328 


0.109 


20.828 


0.589 


6.480 0.182 


RSFI 


0.743 


0.248 


2.187 


0.148 


21.153 


0.758 


6.417 0.182 


RSFII 


0.608 


0.296 


2.219 


0.147 


20.756 


0.774 


6.365 0.206 


SHa 


0.614 


0.262 


2.125 


0.169 


20.138 


0.806 


6.182 0.228 


E+A 


0.526 


0.328 


2.223 


0.211 


18.922 


1.269 


6.067 0.316 



distribution. The mean values of log cr are 2.228 for QST and 2.063 
for E+A, with standard deviations of 0.135 and 0.192, respectively. 
One might argue that the velocity dispersion measurement of E+A 
galaxies may not be a reliable tracer of their masses. The central 
regions of E+As will contain young stellar populations, and ordered 
motion in these central regions could underestimate the velocity 
dispersions. Figure|9](middle and right panels), however, show that 
the distributions of Mz and stellar mass are consistent with that of 
the velocity dispersion of galaxies, confirming that the small log cr 
is part of the nature of E+A. Detailed mass functions of E+As are 
being investigated by Inami et al. (in preparation). 

The lower right p anel of Figure [8] shows the Kormendy rela- 
tion l Kormendvlll977h . a projection of the FP. Black circles and red 
circles denote QST and E+A, respectively, of the same log a range 
(2.2 ^ log cr < 2.4). Solid lines are the best fits of the sample. E+A 
follows its own scaling relation. The tendency towards a smaller 
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Figure 9. Distributions of velocity dispersion (left panel), z-band absolute mag (middle panel), and stellar mass (right panel) for volume-limited early-type 
subsamples. The volume-limited sample is constructed based on the following criteria; 0.025 ^ ^ ^ 0.09 and Mz -20.5. Their stellar masses are obtained 
from the Garching SDSS catalog ^http://vyvyvt'.mpa-garching.mpg. de/SDSS/DR4/ 1. The black soUd is QST, blue dotted line is RSF II, and red dashed line is 
E-hA. 
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Figure 10. The mean log Re and the mean of the five types of galaxies in each log a bin. We divide cr, 70 km s~^ to 330 km s~^, into four bins with the 
bin size of 65 km s~^: 70 ~ 135, 135 ~ 200, 200 ~ 265, and 265 ~ 330 km s~^. The corresponding logarithm values are approximately 1.85, 2.13, 2.30, 2.43, 
and 2.52. The upper two panels compare QST with RSF II and E-l-A. The lower two panels compare QST with RSF I, SHo and E-hA. The properties of RSF 
II are significantly analogous to those of E-l-A in both comparisons. 
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Figure 11. Mean H<5 EW of QST, RSF II, and E+A in each log a bin. The 
range of log a and the bin size are the same as Fig. llOl Note that RSF II has 
about 1 A larger H5 EW than QST. 



re and the more luminous fj. are apparent: the E+A slope is about 
2 times steeper than that of QST, which we believe is responsible 
for the tilt of the E+A PP. In short, both the central concentration 
of starbursts in E+A and its smaller velocity dispersion play an im- 
portant role in the shift and the tilt of the E+A PP. 

We have compared the E+A PP with the QST, RSP V II, 
and SHq PPs in terms of slope, position, and scatter. We wonder 
whether the four types of early-type galaxies and E+A are inde- 
pendent or instead associated in some ways. The origin of the tilt 
and shift of the E+A PP is also explored by analyzing its struc- 
tural and dynamical parameters. When investigating the relation- 
ship among them, the geometry of SP within E+As (i.e. the degree 
of central concentration) is believed to be an important key. We 
discuss the possible evolutionary connections among the different 
types of galaxies in the following section. 



4 DISCUSSION 

4.1 Possible Connection between E+A and UV-excess 
Galaxies 

The PP analysis can provide important clues to the physical links 
among the five different types of galaxies. We have found that there 
are systematic differences in position, slopes and scatter among 
them (see Section 13. 2t . We have shown that the coefficients in- 
creased systematically from SHq (0.55), to E+A (0.63), to RSP 
I (1.03) & II (1.06), and finally to QST (1.12) in the g-band. The 
systematic change may imply a relationship among the different 
types of galaxies. Our main goal here is to examine whether there 
are any evolutionary paths from SHq to RSP I, RSP II, or QST via 
E+A. 

We explore the distribution of the PP parameters (see TableO. 
The largest mean effective radius found in RSP I, despite its rela- 
tively small mean velocity dispersion, suggests that weak ongoing 
star formation in RSP I is not centrally concentrated, but rather 



widely spread across the galaxy. The mean surface brightness of 
RSP I is fainter than that of QST, which also supports the pos- 
sibility. RSP II, however, has a smaller effective radius and more 
luminous surface brightness than QST, which may lead the RSP II 
PP to be analogous to that of E+A. Contrary to our expectation that 
SHq is the most similar to E+A in structural and dynamical param- 
eters, SHq has a significantly larger mean effective radius and less 
luminous mean surface brightness than E+A. These results indicate 
that the ongoing starburst in SHq is not as highly concentrated as 
in E+A (~ 75 per cent of E+As show a positive colour gradient 
according to Yamauchi & Goto 2005), and so the two may have 
significantly different RSP geometries. A sudden quench of RSP 
in SHq would result in E+A with negative or flat colour gradients, 
rather than E+A with blue cores or positive colour gradients. 

PigurefTOlcompares the mean surface brightness and the mean 
effective radius of each galaxy type in the same log a bins. We 
break a, 70 km s^^ to 330 km s^^, into four bins of size 65 km s^^: 
70 ~135, 135 ~ 200, 200 ~ 265, and 265 ~ 330 km s-\ They cor- 
respond approximately to log a = 1.85, 2.13, 2.30, 2.43, and 2.52. 
The bin size of 65 km s~^ is determined in order to generate a 
fair number of log a bins while retaining statistically meaningful 
number of galaxies in each bin. Por clarity, we plot QST, RSP II, 
and E+A in the upper panels and QST, RSF I, SHq, and E+A in 
the lower panels. The short vertical error bars denote the estimated 
uncertainties based on the bootstrap method. In left panels in Pig- 
ure[TO] we compare the mean surface brightness. RSP II, E+A, and 
SHq are more luminous than QST. RSP I, however, shows less lu- 
minous mean surface brightness than QST in the first three bins. 
The same trend is observed for the mean surface brightness of RSP 
II and E+A: an increase in the first three bins followed by a decrease 
in the last bin. Unlike RSP II and E+A, the mean surface brightness 
of SHq remains roughly constant with log a in the first three bins 
and increases in the last bin. The right panels in Pig. [TO] show that 
RSP II and E+A have a smaller mean effective radius than QST in 
all bins, whereas the mean effective radius of RSP I and even SHq 
are larger than that of QST in some bins. 

We found no evidence for connection of RSP I or SHq to E+A 
in this analysis. It is now clear that even though SHq is found in a 
very similar region to the E+A PP, there is no evidence supporting 
their connection. Although both RSP I and SHq show signs of on- 
going star formation, it is hard to link them without exact informa- 
tion regarding star formation histories. The residual star formation 
in the outskirt of a galaxy can persist after the extinction of central 
starburst (due to the gas exhaustion and/or the AGN feedback). In 
this process, SHq could turn into RSP I. Or, it might be that RSP I 
galaxies are ones that somehow undergo star formation across the 
entire galaxy. In this case, RSP I may or may not be related to SHq. 
In order to understand what causes the RSP I population, we plan 
to investigate SHq and ultra-luminous infrared galaxies (ULIRGs) 
in the forthcoming paper (Bae et al., in preparation). 

RSP II displays similar trends to E+A in both mean surface 
brightness and mean effective radius, and the two structural param- 
eters of RSP II fall between those of E+A and QST. These results 
indicate that the star formation geometry of RSP II is also concen- 
trated in the central region as often observed in E+A, even if the 
influence of the young stellar population is weaker compared to 
E+A. A detailed investigation of the physical geometry of RSP ac- 
tivity and its impact on the Scaling Relations will be presented in 
Choi et al. (in preparation). 

Comparisons of the PPs and their parameters for the five types 
of galaxies suggest that RSP II may have a more intimate link to 
E+A than other types. The Balmer absorption strengths can reveal 
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Figure 12. The g,T,i composite SDSS images for QST galaxies in four different log cr bins. Tiie images are sorted by the velocity dispersion: the 1st row 
(70 ~ 135 km s~^), the 2nd row (135 ~ 200 km s^^), the 3rd row (200 ~ 265 km s~^), and the last row (265 ~ 330 km s~^). Selected are galaxies with 
redshift ranging between 0.05 and 0.10. In each row, the 1st (last) column has lowest (highest) redshift. 



a possible transition between E+A and QST via RSF II. FigurelU 
presents the mean H(5 EW at each log a bin for RSF II, QST, and 
E+A. The mean uncertainty of each type (vertical error bars), es- 
timated by the bootstrap method, is as small as the symbol size. 
As inferred from Figure [H H5 EW of RSF II is in between those 
of QST and E+A, supporting possible links between RSF II and 
E+A. Independently, another useful sanity check on their associa- 
tion comes from visual inspection of the E+A galaxies. If there is 
a connection between E+A and QST classes via the RSF II cate- 
gory, then there should be morphological similarities which might 
not be picked up by the FP parameters. Figures [T2lll3l andll4lcom- 
pare the SDSS images of each class in four different sigma bins. It 
is apparent from the figure that the E+A galaxies do look bulge- 
dominated and show similarities to QST and RSF II classes. This 
is also consistent with RSF II being linked to E+A. 

Figure [15] shows possible evolutionary paths from E+A 
to QST via RSF II by using a stellar population model 



( iBruzual & Charlotll2003h . We combine a representative old pop- 
ulation (the median values of observed H(5 and NUV — r of QST; 
large cross) with a young population of solar metallicity {Z = 0.02). 
The stellar population model grids with added young population 
fractions of / = 0.1, 1, 3, and 10 per cent are plotted for the differ- 
ent ages of t = 0. 1, 0.2, 0.5, and 1 Gyr at fixed metallicity of Z = 
0.02, as indicated by the labels. The stellar population models pro- 
vide two reasonable explanations for the relationship among QST, 
RSF II, and E+A. If the star formation activities were quenched 
quite a while ago, RSF II can be considered as the results of pas- 
sively evolved E+A ("the old E+As"). However, if star formation 
was weak, RSF II represents a minor version of E+A ("the mild 
E+As"). We here refer to RSF II as "E+a", named after "E+A". The 
new acronym represents elliptical galaxies ("E") with a minority of 
A-type young stars ("a"). 

Comparison of observational data with stellar population 
models suggests two possible ways of connections between RSF 
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II and E+A. Although the origin of RSF II can be explained by 
both the "old" E+A scenario and the "mild" E+A scenario, E+A 
has a n important p roperty that should not be overlooked. Accord- 
ing to|Got3 ([ioO^), E+As are very rare objects in the lower-redshift 
Universe. In the "old" E+A scenario, even though E+As are likely 
to evolve into RSF II with time, there are not enough E+As to sup- 
ply the entire RSF II population at any given point in time. The 
scenario predicts that the mass distributions of E+A, RSF II, and 
QST should be similar because the difference among these species 
lies only on the recent star formation rate. According to Figure (9] 
however, E+A, RSF II, and QST have a systematic difference in 
there masses and luminosities. These indicate that a large fraction 
of RSF II galaxies should be explained by the "mild" E+A scenario. 

It is important to estimate the number fractions of RSF II 
galaxies following two different paths. In Figure[T5] the stellar pop- 
ulation model suggests that RSF II galaxies right below QST group 
(i.e. -3 A ^ H(5 EW < 2 A) fit better into the "mild" E+A scenario. 
On the other hand, the rest of RSF II galaxies located at the plume- 
like region toward E+As are better explained by the "old" E+A sce- 



nario. Defining the old E+A regime as 2 A (the right edge of the 
mild E+A group) ^ H5 EW ^ 4 A (the left edge of the E+A re- 
gion), we estimate the number of RSF II galaxies belonging to this 
group to be 103 (~ 5 per cent) among 1,936 RSF II galaxies. We 
also estimate the number fraction of E+As in our sample. For a fair 
estimate we only consider the number of E+As with GALEX UV 
information and the number fractions of QST, RSF II, and E+As 
turn out to be 1.0 : 0.9 : 0.04. The number of E+As (~ 5 per cent) 
is comparable to that of RSF II of the plume-like region, further 
supporting their "old" E+A origin. We therefore conclude that a 
large fraction of RSF II (~ 95 per cent) should be the minor (i.e. 
weaker star formation) versions of E+A (i.e. mild E+A), rather than 
passively evolved E+A (i.e. old E+A). 



5 SUMMARY AND CONCLUSION 

The E+A galaxies may undergo the transition from 'blue cloud' to 
'red sequence' and eventually migrate to red sequence early-type 
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galaxies. An observational validation of this scenario is to identify 
the intervening galaxy population that smoothly bridges the E+As 
to red-sequence early-type galaxies. In order to investigate the pos- 
sible associations between E-l-As and UV-excess galaxies, we have 
compared their Fundamental Planes. We used the largest sample 
of 1,021 E-l-As selected from the SDSS DR6 and ~ 20,000 early- 
type galaxies with GALEX UV data. Based on recent star forma- 
tion modes derived from NUV — r colour and Ha emission line, 
we have broken early-type galaxies into QST (RSF-free), RSF I 
(ongoing weak SF), RSF II (post weak star formation), and SHq 
(ongoing starburst) types. The E+A FP is different from the QST 
FP, and this is most likely due to a recent starburst in the central 
region and their intrinsically smaller velocity dispersion. Our sam- 
ple indicates that RSF II galaxies underwent starbursts that were 
weaker than those observed in E-l-As ("the mild E-l-A scenario") or 
the products of passively evolved E-l-As ("the old E+A scenario"). 
The galaxies are characterized by UV-excess but no Ha emission, 
and this is essentially a conceptual generalisation of "E-l-A", in that 
the Balmer absorption line in the "E-l-A" definition is replaced with 



UV - optical colours that are far more sensitive to RSF than the 
Balmer lines. We refer to these UV-excess galaxies as "E-l-a" galax- 
ies (named after "E-l-A"), which stands for elliptical ("E") galaxies 
with a minority of A-type ("a") young stars. We suggest that most 
old "E-l-A" galaxies represent the most recent arrivals to the red se- 
quence. 



A comparative study between E-l-a and E-l-A in terms of their 
preferred environment will give additional insight into their rela- 
tionship. For instance, E-l-As are preferentially found in low-density 
environment in the n earby Univer se and have an excess of neigh- 
bours on small scale ( lGotdl2005bl) . A prediction from the present 
study is that the E-l-a population would also share these properties 
with E-l-As. Schawinski et al. ( 2007) showed that the RSF fraction 
is higher in the field environment, although they do not discrimi- 
nate between RSF I and RSF II. Further investigation is in progress 
on the environment of E-l-a and E+A galaxies in terms of their local 
galaxy number density (Choi et al., in preparation). 
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Figure 15. H5 vs. NUV — r for QST. RSF II, and E+A, overlaid with the simple stellar population model grid. A median SED of QST is used to represent 
the old population. We then add the young population (Z = 0.02) with various combinations of age and mass fraction to the underlying old population. 
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